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We propose a novel generalized method for mass measurements based on phase space singularity 
structures that can be applied to any event topology with missing energy. Our method subsumes the 
well-known end-point and transverse mass methods and yields new techniques for studying "missing 
particle" events, such as the double chain production of stable neutral particles at the LHC. 
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Introduction . At the Large Hadron Collider (LHC), 
physics beyond the Standard Model (SM) may reveal it- 
self in signals with large missing transverse energy. The 
reason is that many TeV scale models require new sym- 
metries that distinguish the new particles from the SM 
fields, which can result in stable neutral particles that can 
be attractive dark matter candidates. One example is the 
neutralino lightest supersymmetric particle (LSP) in the 
minimal supersymmetric Standard Model (MSSM). The 
cascade decays of superpartners result in final state LSP's 
that are "missing particles" that escape detection. 

Measuring the masses of these missing particles is of 
great importance but is challenging. The full S-matrix 
can be used to determine the masses through the event 
distribution profile, but this is highly model-dependent. 
Mass measurement techniques thus rely on kinematic 
analysis, which does not allow for global fits using the 
event profile, but provide useful information via the phase 
space structure defined by the kinematic constraints. 

Event topologies with missing energy may or may not 
be directly reconstructable. For reconstruct able events, 
the number of constraints equals or exceeds the number 
of invisible particle momentum components. By recon- 
structing the momenta, the likelihood of a given test mass 
parameter for each individual event can be obtained 
Though powerful, this method can be used only for exclu- 
sive processes and there are typically many combinatoric 
factors due to the large number of particles. 

For non-reconstructable processes, it is not possible to 
obtain a likelihood contribution from individual events or 
fit the event profile globally using only kinematics. Mass 
measurement techniques include using the end-points of 
kinematic variables such as the invariant mass distribu- 
tion [2], or kinematic cusps [3]- Recently, there has been 
an emphasis on "implicit" variables that depend on trial 
masses [4| , such as the end-point of the mT2 distribution, 
which has a kink when the trial masses equal the true 
masses [5] . Momentum reconstruction is also possible for 
events near the mT2 end-point [5] , and there are attempts 
to understand the mT2 kink based on end-points [7]. 

It is not an accident that these methods use end-points, 
cusps, and kinks, which are singularities in the observable 
phase space. In this letter, we develop the general the- 
ory of kinematic singularities and provide a systematic 



method for obtaining new implicit optimized variables 
that best exploit the singularity structure. We begin with 
a systematic analysis of phase space singularities. We 
then construct kinematic variables that we call singular- 
ity coordinates and apply the method to cascade decays 
and double missing particle chains to show how previ- 
ously studied examples are unified within our approach. 
The mathematical details and additional examples will 
be presented in a future publication [8]. 

Kinematic Singularities. A singularity is a point where 
the local tangent space cannot be defined as a plane, or 
has a different dimension than the tangent spaces at non- 
singular points. The full phase space does not exhibit sin- 
gularities (we do not consider cases of singularities due 
to soft or coUinear massless particles). However, if only a 
subset of the momenta are measured, the relevant quan- 
tity is the projection of the full phase space on to the 
observable phase space of the measured momenta. Each 
observable phase space point can correspond to multi- 
ple configurations in the full phase space, as in Fig. [T] 
The multiplicity volume changes abruptly upon crossing 
a boundary where phase space folding occurs. The ob- 
servable phase space then has a non-smooth density, rem- 
iniscent of caustics in optics. 

Given the visible momenta qj{j = 1, . . . ,rn) and the 
invisible momenta Xk{k = 1, . . . ,n), the full phase space 
is embedded in (n -I- m)-dimensional Euclidean space as 
the solution space of the N constraints: 

g,{x,q) = 0, {i = l,...,N), (1) 

where the gi are coupled polynomial equations that are at 
most of quadratic degree. In mathematics terminology, 
such a space is called an affine variety. 

At the singularity, at least one direction of the tangent 
plane in the full phase space is aligned vertically along the 
invisible momentum direction. The invisible momentum 
components of the normal vector that defines this tangent 
plane are given by the row vectors of the "restricted" 
Jacobian matrix (dgi/dxk). The vertical alignment of the 
tangent plane implies that at a singularity, the restricted 
Jacobian matrix must have a reduced rank: 

Rank ( ^) < Rank ( ^] . (2) 
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FIG. 1; A schematic diagram describing the relation between 
the full phase space and the projected observable phase space. 

The amount that the rank is reduced is the degree of the 
singularity; a wall (cusp) has degree one (two) . 

Finding the reduced rank condition of an arbitrary ma- 
trix is not an easy problem. However, for the special case 
of an affine variety, the given set of polynomial equations 
can be substituted by another set with the same solution 
space. The set of all such equivalent equations is called 
an ideal, which is generated by a finite set of polynomials 
known as a basis. We focus here on the Grobner basis, in 
which variables are eliminated sequentially as follows: 

gi{xi,X2,X3, . . . ,x„) = 0, 

g2iX2,X3, ...,Xn) = 0, 
gNixN,XN+l, ■ ■ ■ ,Xn) = 0. (3) 

The algorithm for finding the Grobner basis for a general 
coupled polynomial system is known P]. For the pro- 
cesses of interest, it is tractable to obtain it analytically. 

The reduced rank condition implies that one or more 
row vectors of the restricted Jacobian are linearly depen- 
dent. In the Grobner basis, the restricted Jacobian is of 
upper triangular form. Therefore, a necessary but not 
sufficient condition for linear dependency is that one of 
the diagonal components vanishes, resulting in an ana- 
lytic condition for the singularity position. 

Singularity Coordinates. The next step is to construct 
an optimized one-dimensional variable that we call the 
singularity coordinate. This is an implicit variable be- 
cause the location of the singularity is defined by the re- 
duced rank condition of the restricted Jacobian matrix, 
which is an implicit function of the mass parameters. The 
singularity coordinate must satisfy the following criteria: 
(i) it must be zero at the singularity, (ii) its direction 



must be perpendicular to the singularity hypersurface in 
observable phase space, and (iii) it must be normalized 
such that every event can give the same significance. 

To see this, we note that the reduced rank condition 
implies that one linear combination '^^Ci{dgi/ dxk) be- 
comes a null vector at the singularity point. The perpen- 
dicular direction is determined by {v)^ = Ci{dgi/dqj), 
(recall qj are the visible momenta) . The singularity coor- 
dinate in this direction is maximally efficient for reveal- 
ing the singularity structure. To assign an unambiguous 
value to each event, the singularity coordinate is scaled 
so that events with the same invisible phase space volume 
around the nearest singularity have the same value. This 
requires a knowledge of the local phase space properties 
around the singularity at quadratic order. 

As shown in Fig. [2j a local orthonormal coordinate 
system around a given reference point can be split into 
tangent directions tr,{r — 1, ... n 4- to — N) and nor- 
mal directions n^, (s = 1, . . . , N). A general phase space 
point near this reference point is labeled by the tangent 
coordinate. The normal coordinate is determined by a 
quadratic function of the tangent coordinate (the second 
fundamental form): = IP{t^) = M^^it^t^ , where 

We define S = t/ • II{tr) = Mrr'trtr' . To find the ap- 
propriate scale factor, we need to obtain the phase space 
volume in the tangent directions that correspond to in- 
visible momenta. The phase space in the invisible mo- 
mentum directions in the diagonalized basis t^ is given by 
01^1+' ■ ■+flAf^M — where M is the number of invisible 
tangent directions and the eigenvalues determine the 
shape of the invisible phase space around the singular- 
ity. For positive eigenvalues, the ellipsoid-shaped phase 
space volume scales as (Vol) oc (aia2 . . . OAf )"^^^ S^^^^. 
The singularity coordinate S that satisfies all three cri- 
teria is thus given by 

See {ai . . .omT^'^'^ ^- (5) 

S is an implicit kinematic variable, since the location of 
the zero, the normal direction v and the scale factor can 
be defined only when mass parameters are given. 

Simple Cascade Decay. Our first example is the sim- 
ple cascade decay process shown in Fig.|3[(e.g. neutralino 
decay in the MSSM, with Y = xl, L = 1, and X = Xi). 
The on-shell equations of this system are 

x^^rux, [x + qff = m\, {x + qj + q^Y = m\ . (6) 

"nix, 'niL, and tuy are trial masses, x is the invisible parti- 
cle momentum, and qnj are the visible particle momenta. 
Taking the z-axis in the direction of the 3-momentum of 
In in the center of mass (CM) frame of the visible par- 
ticles In.}, we have = (£'cm/2, 0, 0, ±i5cin/2), where 



FIG. 2: The scaling behavior near a singularity. 
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FIG. 4: The event topology of a double missing particle chain. 





FIG. 3: The event topology of a simple cascade decay. 



The scale factor can be read off from 53 = 0. The sin- 
gularity coordinate is given by the well-knovifn quantity 
S cx ((mJ7^'^)2 - £:2^)/(m^^'^)2. Our method also shows 
that no other singular structures exist for this process. 

Double Missing Particle Chain. For the double miss- 
ing particle chain of Fig. |4j we have the relations 



Ecm is the CM energy of Inj- Eq. (|6| is then 



(£;„„/2 + xof -x\-xl- + x^f 

(E'cm + a^o)^ 



x\ 



m,Y- 



(7) 



The Grobner basis for this system is particularly simple. 
With the lexicographic ordering xq X3 Xi X2, 
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The restricted Jacobian (dgi/dxj) is 



2Er 



2El^^Xi 2El^X2 



(8) 



(9) 



The first two row vectors are zero only when there are soft 
singularities, which we do not consider here. The condi- 
tion that the third vector vanishes results in xi — X2 — 0. 
Physically, this means that the missing particle momen- 
tum is aligned in the direction of the lepton momentum 
in the CM frame. Together with Eq. ([t]), this results in 
the following condition at the singularity: 



cm o - 



{m^r^)\ (10) 



{x2 + q2f = m 



2 / , \2 2 

ix, (Xi+qi) =mY, 

Y, XiT + X2T=PT: 



(11) 



{xio, Xii, Xi2, Xi^) denote the momenta of Xi, 
= (giOj fti: Qi2,<li3) are the visible particle momenta for 
each chain, and pt = {pti,PT2) is the missing transverse 
momentum. We assume mxi — mx2 — "nix and myi — 
mY2 = my (the generalization to an asymmetric chain is 
straightforward). With xio >- 2^13 >- X20 >- X21 >~ X22 >- 
X23 >~ Xii >- X12, the Grobner basis takes the form 
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= gioa^io - 9132:13 - giia;ii - gi2a;i2 - Ci, 



(12) 



(910 - 913) 2^13 - 2giigi3Xi3a;ii - 2gi2(7i3a;i3a;i2 



-2Cigi3Xi3 -I- (qIq 



9?i) 



"11 



29ii9i2a;iia;i2 





+ (910 - 9^2) ^12 


- 2CigiiXii - 2Cigi22:i2 




+ {mWi. - Cl) 


(13) 
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= 9202:20 - 9232:23 4 


- 9212:11 + 5222:12 - (^2, (14) 
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= 2:21 + .T11 -PTI, 


(15) 
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= 2:22 + .T12 -PT2, 


(16) 
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(2 2 \ 2 . 
— 1920 ~ 923j ^23 + 


29219232:232:11 + 29229232^232:12 



-2(729232:23 + (920 - 92l) 



^11 



29219222:112:12 



+ (920 - 922) XI2 - (2pTi920 - 2C2921) 2:11 (17) 



-{2p 



'T2920 



2C2922) 2:12 + {Pt + ITT-x) 920 ~ Q 



2 J 



This reproduces the well-known result for the edge of the 
invariant mass mu = \/ (9n + 9/)^- The tangent direc- 
tions here are given by 2:1^2, and v is the Ecm direction. 



in which Ci — (my — m^ — 91 • 9i)/2, C2 = {rriY — m,3 
92'92+92T'Pt- The restricted Jacobian matrix (dgi/dxj) 
with respect to the invisible momenta x^ 2 has the form 



4 



4 



^ 600 
!?500 



I 400 



300 



200 



100 



m'™^ = 200 GeV, r^™^ = 500 GeV 
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FIG. 5: The singularity coordinate distribution for several 
choices of trial masses with mx ~ 200 GeV, my = 500 GeV. 
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in which □ is a nonzero entry that depends on the visible 
and invisible momenta, 1 is a constant (nonzero) term, 
and a blank space is a zero. The conditions dgi/dxiQ — 
qiQ and dg3/dx2o — 920 correspond to soft singularities. 
The restricted Jacobian (dgi/dxj) has a reduced rank if 



dg2 
dxi3 



dx23 



2(<7io ~ 913)2^13 - 2(Ci + qiT ■ xiT)qi3 = 0, 



2(920 - 923)2^23 - 2(C2 - 92T • ^lT)g23 = 0, 
892 dg-2 



(18) 



Once the nearest singularity point is identified, we can 
determine the singularity coordinate. A numerical anal- 
ysis for events with true masses mx — 200 GeV and 
my = 500 GeV is shown in Fig. [5] The singularity ap- 
pears at S = only when the trial masses are equal to 
the true masses, providing a proof of concept. 

Here we have neglected backgrounds and assumed that 
we have identified the event topology correctly. A com- 
plete treatment of the backgrounds must be done on a 
case-by-case basis and is beyond the scope of this paper. 
However, as the singularity coordinates maximize any 
singular features at the true masses, they best discrim- 
inate the signal in the presence of backgrounds, which 



have smooth profiles at these points. If the event topol- 
ogy was misidentified, singular features will not appear 
for any trial mass values, indicating that the hypothesis 
was incorrect. Further details will be given in [5]. 

Conclusions . We have presented a systematic method 
for measuring missing particle masses based on phase 
space singularity structures that is applicable to any non- 
reconstructable process. The approach reproduces well- 
known results, such as the end-point of the invariant 
mass distribution, and provides a means for finding the 
singularities for more general processes that can also be 
used for determining qualitative event properties. This 
method should provide invaluable tools in the search for 
new physics at the LHC. 
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